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SU(3) vortex- like configurations in the maximal center gauge 
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A new algorithm for fixing the gauge to (direct) maximal center gauge in SU(N) lattice gauge theory is presented. 
We check how this method works on SU(3) configurations which are vortex-like, and show how these configurations 
look like when center projected. 
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1. Introduction. 

The maximal center gauge is widely used to 
study low energy phenomena such as the confine- 
ment property or the breaking of the chiral sym- 
metry, as can be seen in a set of recent works 
jl]-^ and also in these proceedings. Nevertheless, 
at least to our knowlegde, there was no efficient 
method of direct gauge fixing to maximal center 
gauge in SU(N) lattice gauge theory for values 
N > 2. This is the motivation of the first part 
of this work in which we present a new algorithm 
of direct center gauge fixing to maximal center 
gauge. We apply this new method to previously 
prepared SU(3) vortex-like configurations. Our 
purpose in the second part of this work is to know 
how these configurations look like in the maximal 
center gauge, and therefore, whether these solu- 
tions have the properties described in references 
jl]-^ for a confining object. For a more detailed 
description of our results see reference . 

2. The method. 

The maximal center gauge in SU(N) lattice 
gauge theory is defined as the gauge which brings 
fink variables U as close as possible to elements of 



the center Zn = {e^'^™*/^! 
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This can be achieved by maximizing the following 
quantity: 

R= ^ I Tr C/„(n) |2 (1) 

where NsUe is the number of sites in the lattice 
and Ndim the number of dimensions {R satis- 
fies \R\ < 1). Our procedure is based on a lo- 
cal update of R. We maximize this quantity re- 
spect to a gauge transformation G(n) defined at 
the lattice point n. We build the SU(N) ma- 



trix G(n) from a SU(2) matrix g(n) as in the 
Cabbibo-Marinari-Okawa method. Doing the up- 
date in this way the problem is reduced to find- 
ing the maximum of a quadratic form, and this 
can be done using standard methods (see (lOj 
for details). Once we obtain the SU(2) matrix 
g{n) maximizing i?, we update the t/^ variables 
touching the site n: U^{n) G{n)Uf_i{n) and 
U^{n — fx) — > Ufj,{n — fi)G^n). We repeat this 
procedure over the N{N — l)/2 SU(2) subgroups 
of SU(N) and over all lattice points. When the 
whole lattice is covered once we say we have per- 
formed one center gauge fixing sweep. We make 
a number of center gauge fixing sweeps on a lat- 
tice configuration and stop the procedure when 
the quantity R is stable within a given precision 
(10~^ in this work). 
3. The vortex solution. 

Following the same methods used in |0 for the 
SU(2) group, we built a SU(3) Yang-Mills config- 
uration which has the following features: 

• It is a solution of the SU(3) Yang-Mills clas- 
sical equations in i?^. 

• It is constructed from R^ x solutions of 
the SU(3) Yang-Mills classical equations by 
glueing to themselves the two periodic di- 
rections {R^ X is considered the limit of 

with two periods much larger than the 
other two). 

• The R^ X building block satisfies non- 
orthogonal twisted boundary conditions, is 
(anti)self-dual and has minimal action S = 
87r^/3. The crucial point for this building 
block to have vortex-like properties is that 
the twist in the small plane must be non 
trivial. 
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We isolate the x building block us- 
ing the lattice aproach. Working on lattices of 
sizes Ni X Ni X Ns X Ng (directions t,x,y, z re- 
spectively and Ni» Ng) and imposing twisted 
boundary conditions given by the twist vectors 
fc = m = (1, 0, 0) we isolate SU(3) configurations 
with topological charge Q = —1/3 and minimal 
action S = Sn^/S using a standard cooling algo- 
rithm. The lattice sizes used are Ng = 4, 5, 6 and 
Ni = 6 X Ng, and we obtain configurations with 
action 5/(87rV3) = 0.92219,0.95015, 0.96536 re- 
spectively. We fix the length in the small direc- 
tions equal to 1 and then the lattice spacing is 
a^l/Ng. 

4. Properties of the solution 

A) From the field strength F^i,, obtained using 
the clover average of 1 x 1 plaquettes, we calcu- 
late the quantity £^iy{x^, x,^), the integral of the 
action density in two directions, defined as, 

£t,,^ixfj.,x„) = j Sxa^^^u ^Tr (Fpo-Fp,^) (2) 

and obtain the following results: 

• Localized in the xt plane as shown in the 
plot of £xt{x,t) included in figure 1. 

• Localized in the x direction and almost flat 
in the y direction as shown in the plot of 
£xy{x,y) included in figure 1. In this case 
we repeat six times the y direction to obtain 
a box size equal to the one used in the figure 
for £xt{x,t). Identical figures are obtained 
changing x ^ t and y ^ z. 

• Almost flat in the y, z directions as shown in 
the plot of Syziy, z) included in figure 1. We 
repeat six times both directions to obtain a 
box size equal to the one used for Sxtix, t). 

B) We calculate the Wilson loop Wc{r) around 
this object, 

Wc{r) = 1 Tr (^Pexp tA^dXf}j (3) 

being C an r x r square loop in the xt plane 
centered at the maximum of the solution. We 
parametrize Wc{r) by the functions L{r) (its 
module) and 0(r) (its phase). In figure 2 we 
show the quantities L{r) and 0(r) as a function 
of r, putting in the same plots the values with the 
square loop centered at the minimum of the solu- 
tions in the z directions and at the maximum 




Figure 1. Action densities integrated in two di- 
rections. 

in the t, x directions. The conclusions extracted 
from figure 2 are the following: 

• These functions arc almost independent of 
coordinates y and z. 

• When r is bigger than the size of the object 
we obtain Wc{r) ~ ea;p(i47r/3), an element 
of the group center. 

These are the two expected properties for a vor- 
tex. 

5. The center-projected configuration. 

We apply to the obtained solutions the method 
of center gauge fixing. For the SU(3) vortex- like 
solutions we have worked with we have obtained 
the following maximum values: R = 0.90123970, 
0.93246431 and R = 0.94918155, for the lattices 
sizes: iVg = 4, 5 and Ng — 6 respectively. We 
project the SU(N) link variables to Z^v link vari- 
ables and we calculate the values of the plaque- 
ttes from the link variables. We obtain the 
following results: 

• In the xt plane, the plane in which the vor- 
tex is localized, we obtain for the configu- 
rations with maximum value of R the same 
structure for all y, z points. Only two pla- 
quettes are different from the identity. The 
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Figure 2. The module L{r) (top figure) and the 
phase (j){r) (bottom figure) of the Wilson loop de- 
fined in equation]^ as a function of r. Full symbols 
correspond to Wilson loops centered at the max- 
imum of the solutions and empty symbols to the 
same quantity centered at the minimum of the so- 
lutions in the y, z directions and at the maximum 
in the t, x directions. 

first one is located at the top-right corner 
and reflects the use of twisted boundary 
conditions. The other one is located near 
the maximum in the action density of the 
solution and has the same value of the Wil- 
son loop shown before when the size of the 
loop is much bigger than the size of the so- 
lution. This is one of the most interesting 
results of our work, showing how the vor- 
tex properties are reflected in the center- 
projected configuration. 

• In the xy plane ( or xz,ty, tz plane). In this 



plane the vortex is localized in one direc- 
tion and almost flat in the other. There is 
no regular structure for all t, z points. The 
most common configuration is that one with 
all plaquettes equal to the identity. 
• In the yz plane. The solution is almost 
flat in both directions. As in the previous 
case, there is no regular structure for all t, x 
points. The most common configuration is 
that one with a very similar structure to the 
one in the xt plane. 
6. Conclusions. 

The main conclusions of this work are the fol- 
lowing two. First, our method of gauge fixing 
to maximal center gauge works very well. And 
second, we have seen that, for the configurations 
with maximum value of R, there is a clear relation 
between the structure of the SU(3) configuration 
and the structure of the center-projected one. 
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